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ABSTRACT 



> 

CN . Observations of the optical afterglow of GRB 041006 with the Kiso Observa- 

^ . tory 1.05 m Schmidt telescope, the Lulin Observatory 1.0 m telescope and the 

Xinglong Observatory 0.6 m telescope. Three-bands (B, V and R) of photometric 
data points were obtained on 2004 October 6, 0.025—0.329 days after the burst. 
These very early multi band light curves imply the existence of a color dependent 



plateau phase. The 5-band light curve shows a clear plateau at around 0.03 days 



after the burst. The R band light curve shows the hint of a plateau, or a possible 
slope change, at around 0.1 days after the burst. The overall behavior of these 
multi-band light curves may be interpreted in terms of the sum of two separate 



components, one showing a monotonic decay the other exhibiting a rising and a 
falling phase, as described by the standard afterglow model. 
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Subject headings: Gamma-ray Bursts: afterglow 

1. Introduction 

Based on the standard afterglow model (e.g. Sari, Piran and Narayan 1998), an optical 
light curve can be expected to consist of a combination of four power-law segments, with 
connections at certain break frequencies. The break frequencies are the self-absorption fre- 
quency (usa), the typical frequency (u m ) and the cooling frequency (u c ). Before the HETE-2 
era, almost all GRB light curves showed only the v > v c segment characterized by a single 
declining or achromatic break. Note that the standard afterglow model predicts the existence 
of a peak at an early time when the typical synchrotron v m frequency crosses into the optical 
frequency. Before the peak time {y ov t < v m ), the luminosity increases proportional to t 1 ^ 2 , 
until reaching the maximum flux -F max at v — u m ; the subsequent decay of the luminosity 
is proportional to £ 3 ( 1_ p)/ 4 , where p is the index of the power-law distribution of energetic 
electrons accelerated at the shock. Since the break time varies with the frequency, multi- 
color observations of the very early afterglow are required, so as to catch u, m in the optical 
wavelength. Up to now, there have only been a few cases (such as GRB 021004; Urata et al 
2005a) for which early time multi-color data are available. 

GRB 041006 was detected with HETE-2 at 12:18:08 UT on 2004 October 06. The Wide- 
Field X-ray Monitor (WXM; Shirasaki et al 2003) localized the burst in real time, resulting 
in a GCN alert 42 seconds after the burst trigger. The flight error region was a circle with a 
14' radius (90% confidence level) centered at a 2000 = 00 h 54 m 54 s , 5 2000 = 01°18'37". Spectral 
analysis showed the 2-30 keV fluence of this event to be 5 x 10~ 6 erg/cm 2 and the 30-400 
keV fluence to be 7 x 10 _6 erg/cm 2 : the ratio between these two values means that it can be 
classified as an "X-ray rich GRB" . The light curve shape of the gamma-ray pulse of GRB 
041006 is very similar to that of GRB 030329 but, its spectral characteristics are 20 times 
fainter. GRB 041006 shows a soft precursor before the main gamma-ray pulse (Galassi et 
al. 2004). At 1.4 hours after the burst, the optical afterglow was found within the 14'-radius 
error circle at the following coordinates: a 2000 = 00 h 54 m 50 s .17, <5 2000 = 01°14'07" (Da Costa 
et al. 2004). The redshift was determined by Price et al. (2004) using the Gemini North 
telescope to be z = 0.716. 
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2. Observations 

The follow-up observations of the GRB 041006 optical afterglow at the Kiso, Lulin and 
Xinglong Observatories were carried out within the framework of the East Asia Follow-up 
Observation Network (EAFON, Urata et al. 2005b). The 1.05 m Schmidt telescope and a 
2k x 2k CCD camera at the Kiso observatory were used to make 300 sec exposure B, V, and 
R band imaging observations, starting at 12:51 UT on 2004 October 6 (0.009 days after the 
burst). The field of view is 51'. 2 x 51'. 2, and the pixel size is 1".5 square. Before receiving the 
report by Da Costa et al. (2004), we had attempted to cover the entire HETE-2 error region. 
We made B, V, and R band imaging observations of 300 sec exposure each. The afterglow 
was clearly detected in all bands. In order to image simultaneously the optical afterglow and 
several Landolt standard stars in SA92 (Landolt 1992), we pointed the telescope 5 arcmin 
to the east (Figured]). 

The R band observations were performed using the 0.6 m telescope at the Xinglong 
Observatory, starting at 2004 October 6 (14:33 UT; 0.094 days after the burst). The filed of 
view is 10'. 1 x 10'. 1, and the pixel size is 0".47 square. The B and R band observation of the 



afterglow was also performed with the Lulin-One-meter- Telescope (jKinoshita et al. 20051 ) on 
the night of October 6 (0.320 — 0.329 days after the burst). Additional B, V, R and I band 
observations were made by the Maidnak 1.5 m telescope on October 6. The time coverage 
of the observations was from 0.149 to 0.188 days after the burst. 



3. Analysis 

The data reduction was carried out using the standard package NOAO IRAF package. 
We performed the bias-subtraction and flat-fielding correction using the appropriate calibra- 
tion data. We calibrated the flux zero-point of our images through a comparison with the 
SA92 standard star field. We also performed cross calibrations of our photometric results 
using several field stars for which the magnitudes have been calibrated by Henden (2004). 
The difference of photometric zero points between our SA92 and Henden's filed photometry 
is within 0.04 mag. Aperture photometry for all the data was performed using the APPHOT 
package of IRAF. The photometric results are summarized in Table [TJ 
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4. Results 
4.1. Light curves 

The multi-band light curves of the GRB 041006 afterglow are shown in Figure [2J In 
addition to our B, V and R band data, there are several unfiltered observations, the R 
band data points reported in the GCN Circulars (Yost et al. 2004; Ayani &; Yamaoka 2004; 
Fugazza et al. 2004; Monfardini et al. 2004; Misra & Pandey 2004; Kahharov et al. 2004; 
Kinugasa & Torii 2004), and a number of published R and V band data points (Soderberg 
et al. 2005; Stanek et al. 2005) were included in Figure 2. The 5-band light curve showed 
a clear plateau at around 0.03 days after the burst, although its behavior in earlier epochs 
remains unknown. Interestingly, the R band light curve shows a hint of a plateau, or a 
possible slope change, around 0.1 days after the burst. These observational results indicate 
that the B, V and R band photometric points obtained by EAFON play an important role 
in characterizing the temporal evolution of the afterglow. 

As is obvious in the Figured the single power law fit (a = —0.95, x 2 / u = 5.76 with 
z^=56) of the R band light curve is poorly described. It is also unlikely that the V band 
light curve can be described with the single power law (a = —0.72, x 2 l v — 1-54 with u=9). 
On the other hand, the data subsets at t < 0.04 days and t > 0.15 days can be described 
successfully by two separate single power-law models. For the R band, we have obtained 
a = -0.91 ± 0.12 with x 2 /" = 0.30 for v = 2 before 0.04 days, and a = -1.12 ± 0.02 with 
X 2 /y = 1.44 for v = 48 after 0.15 days; for the V band a = -0.59 ± 0.03 with x 2 /" = 0.41 
for v = 3 before 0.04 day, and a = —1.41 ± 0.30 with x 2 l v = 0-19 for v = 1 after 0.15 days 
From this point of view, the overall behavior of these multi band light curves may then be 
understood as the sum of two separate components, one showing a monotonic decay with 
the other having a rising and a falling phase. 

4.2. Color change of the afterglow 

The photometric results were corrected for Galactic reddening, using the reddening 
map of Schlegel, Finkbeiner, & Davis (1998). The Galactic reddening toward the burst is 
E(B - V)= 0.022, which implies a Galactic extinction of A B = 0.094 and A R = 0.058. 
Although the light curves in Figure [2] suggest a significant color change within < 0.1 days, 
our data coverage is insufficient to make a conclusion about such a possibility. Accordingly, 
for simplicity's sake we assume that at t > 0.15 days the power-law spectrum of the afterglow 
has a constant index. The observed color, B — R = 0.83 ±0.28 mag, then indicates a spectral 
index of (3 ~ — 1.24±0.42 in terms of f(u) oc . This is close to the value of j3 = — 1.07±0.02, 
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which is predicted by the spherically symmetric model of Sari, Piran, & Halpern (1999) in 
the regime of v c < v opt in combination with the measured a = —1.118 ± 0.016. 



5. Discussion 

We observed the GRB 041006 optical afterglow from the very early phase (~ 0.03 days 
after the burst) i n three bands. Alth ough the relatively early phase of a dozen afterglows such 



as GRB 990123 (lAkerlof et al. 19991 ) . multi-color observations of very early optical afterglow 



are still rare. In the GRB 021004 case, the very early optical after glow shows a clear a chro 



matic re-brightening phase which peaks was at around ~ 0.07 days (lUrata et al. 2005bl ). On 
the other hand, the afterglow behavior of GRB 041006 depends on its color. The B band 
light curve shows a clear plateau at around 0.03 days after the burst, although its behavior 
in earlier epochs remains unknown. The R band light curve shows the hint of a plateau, or 
a possible slope change, around 0.1 days after the burst. 

Several models, based on the standard afterglow model, can be used to explain the 
variability in the light curve, such as the variable external density (e.g. Lazati et al. 2002), 
refreshed shock (Rees Sz Meszaros 1998; Kumar & Piran 2000a; Sari & Meszaros 2000), 
and patchy shell models (Kumar & Piran 2000b). To explain the present multi-band light 
curves, these models require similar temporal behavior in each band. If the light curve 
of GRB 041006 does not show a color change, it should have two bumps with almost the 
same amplitudes. The patchy shell model is not suitable to explain current light curves in 
this case. This is because the model predicts that the amplitude of the fluctuations in the 
afterglow light cure is expected to decrease with time proportional to the Lorentz factor. 

The suggested second component may be identified with the brightening phase of Sari 
et al's standard afterglow model (1998) which predicts the evolution to consist of four phases 
involving a convex-shaped period; the brightening phase can be expressed with a power-law 
index of 0.5. To describe the plateau phase together with the monotonically decreasing 
first component, we use the single power law, plus a smoothly broken power law function 
expressed as 

p* 

Fv{t) = m^+m) *]' (1) 

with a% fixed at 0.5 after Sari et al. (1998). We obtain an acceptable fit for the R band 
light curve with a = -1.07 ± 0.01, a 2 = -1.49 ± 0.03 and t b = 0.143 ± 0.012 {j^jv = 1.19 
with v = 66). We also fit the B and V band light curves with the same function. For better 
constraints, we fix the ot\ = —1.0 for B band and «2 = —1.4 of the V band. Due to the lack 
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of data points for the first components in the B band, and the former components in the 
V band., these vales come from the result of R band fitting. For the B band, we obtained 
a 2 = -1.47 ± 0.07 and t b = 0.0675 ± 0.005 (x 2 /v = 0.47 with u = 7); for the V band, 
a = -1.08 ± 0.02 and t b = 0.073 ± 0.006 {x 2 /v = 0.83 with v = 8). In Figure [21 these best 
fit functions are superposed on the B, V and R band data points. Based on the equation 
(Sari et al. 1998) 



5.7 x 10^ a #& / V /a Hz > 



(2) 



we can estimate the strength parameter of the magnetic field eg and/or the injected elec- 
trons parameter e e as follows; e l ^ 2 e 2 = 4.1 x 10~ 3 . This value is consistent with the result 
of Panaitescu & Kumar (2001, 2002) which was estimated from late-time multi- frequency 
afterglow observations. 

About the first component, the decay index is a = —0.92 and the back-extrapolated 
power-law funct ion can successful ly account for the point at 14.4 min after the burst observed 
by ROTSE-III (jYost et al. 20041 ) . The temporal index is considere d to be directly rela ted 
to the prompt emission, as obse rved in the case of GRB 990123 flAkerlof et al. 19991 ) or 
GRB 030418 (IRykoff et al. 20041 ) . However, the tempor al index is incon sistent with the 
values (—0.5 or —2) predicted by a reverse shock model (IKobayashi 20001 ). Thus the first 
component is not likely to be an optical flash, such as is the case for GRB 990123. Considering 
the above discussions one of possible explanation is that these emissions come from multiple 
jets and/or sub-jets (Nakamura 2000, Yamazaki et al. 2004). In the multiple jet case, 
the former emissions and the latter components are caused by narrower and wider jets, 
respectively. 



6. Conclusion 

Thanks to our EAFON observations, we have obtained early B, V and R band possible 
re-brightening episode for GRB 041006. While there are several other af terglows similar to 
GRB 050319 (Wozniak et al 2005; Huang et al 2006) and GRB 050525a fclotz et al. 2005h 
for which the early optical behaviors may possibly be similar to that of GRB 041006, their 
physical mechanism is not clear, due to the general lack of color /spectrum information on the 
early optical afterglows. The present work has provided some hints for the current shallow 
decaying of the X-ray and optical behaviors. 
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Fig. 1. — An R band image of the GRB 041006 field taken at the Kiso observatory. The 
circles indicate the afterglow and standard stars. North is up; east is left. 
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Fig. 2. — B, V and R ba nd light curves produ ced base d on the Kis o, Lulin and Beijing 
results together with SSO JSoderberg et al 2005b . MMT Jstanek 2005h and several GCNs. 
The dotted lines indicate the best fit model functions described in the text. The dashed 
lines indicate the model components of the best fit function for the R band lightcurve. 
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Table 1: Log of follow-up observations of the afterglow of GRB 041006. 



Date 


Start Time(UT) 


Delay (days) 


Filter 


Exposure (s) 


mag 




Site 


2004-10-06 


13:13:04 


0.040 


B 


300 s 


X 


1 


18.900±0.055 




Kiso 


2004-10-06 


13:20:05 


0.045 


B 


300 s 


X 


1 


18.965±0.058 




Kiso 


2004-10-06 


13:27:05 


0.050 


B 


300 s 


X 


1 


18.890±0.070 




Kiso 


2004-10-06 


15:07:22 


0.119 


B 


300 s 


X 


1 


19.696±0.074 




Kiso 


2004-10-06 


16:09:02 


0.162 


B 


300 s 


X 


1 


19.997±0.074 


Mt. 


Maidanak 


2004-10-06 


16:27:48 


0.175 


B 


300 s 


X 


1 


20.122±0.078 


Mt. 


Maidanak 


2004-10-06 


19:57:08 


0.320 


B 


300 s 


X 


1 


21.004±0.188 




Lulin 


2004-10-06 


14:38:52 


0.099 


V 


300 s 


X 


1 


19.141±0.117 




Kiso 


2004-10-06 


14:52:56 


0.109 


V 


300 s 


X 


1 


19.316±0.145 




Kiso 


2004-10-06 


16:03:40 


0.158 


V 


240 s 


X 


1 


19.564±0.057 


Mt. 


Maidanak 


2004-10-06 


16:14:57 


0.166 


V 


240 s 


X 


1 


19.677±0.078 


Mt. 


Maidanak 


2004-10-06 


16:33:28 


0.179 


V 


300 s 


X 


1 


19.752±0.074 


Mt. 


Maidanak 


2004-10-06 


12:51:20 


0.025 


R 


300 s 


X 


1 


17.799±0.040 




Kiso 


2004-10-06 


12:58:21 


0.030 


R 


300 s 


X 


1 


17.988±0.063 




Kiso 


2004-10-06 


13:05:23 


0.035 


R 


300 s 


X 


1 


18.152±0.080 




Kiso 


2004-10-06 


14:32:51 


0.094 


R 


120 s 


X 


1 


18.97±0.23 




Beijing 


2004-10-06 


15:16:15 


0.124 


R 


120 s 


X 


1 


18.89±0.24 




Beijing 


2004-10-06 


15:49:32 


0.149 


R 


180 s 


X 


1 


19.211±0.067 


Mt. 


Maidanak 


2004-10-06 


15:54:56 


0.152 


R 


180 s 


X 


1 


19.231±0.075 


Mt. 


Maidanak 


2004-10-06 


16:19:59 


0.170 


R 


180 s 


X 


1 


19.229±0.074 


Mt. 


Maidanak 


2004-10-06 


16:39:30 


0.183 


R 


300 s 


X 


1 


19.544±0.076 


Mt. 


Maidanak 


2004-10-06 


16:46:28 


0.189 


R 


300 s 


X 


1 


19.371±0.156 


Mt. 


Maidanak 


2004-10-06 


20:02:54 


0.324 


R 


300 s 


X 


1 


20.349±0.090 




Lulin 


2004-10-06 


20:09:25 


0.329 


R 


300 s 


X 


1 


20.314±0.092 




Lulin 



